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^ We present preliminary absolute branching fraction measurements of 

leptonic Df decays to n'^i'i^ and t'^i't, and of hadronic decays to 

I K~K'^7r'^, K^K^ and ryvr"'". The results are obtained from a large data 

>• sample collected near the T(45') and T(5S') resonances with the Belle 

^ detector at the KEKB asymmetric-energy e^e^ collider. Wc obtain the 

a^ the branching fractions B{D+ n+iy,,) = (0.528 ± 0.028 ± 0.019)% and 

<^ B{Df T+Ur) = (5.70±0.2llg|J)% which are combined to determine the 

^ Dg decay constant fo^ — (255.0 ± 4.2 ± 5.0) MeV, where the first and sec- 

^ ond uncertainties are statistical and systematic, respectively. Branching 

^ — I fractions of hadronic decays are measured to be B{D^ — )■ K~K~^7i~^) = 

> (5.06 ± 0.15 ± 0.19)%, B{Df WK+) = (2.84 ± 0.12 ± 0.08)% and 

B{D+ r]7i+) = (1.79 ± 0.14 ± 0.05)%. 
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1 Introduction 

The leptonic decays of mesons provide access to experimentally clean measurements 
of the meson decay constants or the relevant Cabibbo-Kobayashi-Maskawa matrix 
elements. In the Standard Model (SM) the branching fraction for a leptonic decay of 
a charged pseudoscalar meson, such as Z)j", is given by [H |2]: 



where Md^ is the Dg mass, r^^ is its lifetime, mi is the lepton mass, Vcs is the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element between the Dg constituent quarks c and 
s, and Gp is the Fermi coupling constant. The parameter fa^ is the decay constant, 
and is related to the wave-function overlap of the quark and anti-quark. The leptonic 
decays of pseudoscalar mesons are suppressed by helicity conservation and their decay 
rates are thus proportional to the square of the lepton mass. Leptonic Ds decays into 
electrons are not observable whereas decays to taus are favored over decays to muons 
in spite of the reduced phase-space in the former case. 

If the magnitude of the relevant CKM matrix element is well known from other 
measurements then by measuring the leptonic branching fraction of a pseudoscalar 
meson one can determine the decay constant with high precision. Conversely, if one 
can precisely estimate the decay constant of a pseudoscalar meson, it is possible to 
determine the magnitude of the relevant CKM element. 

Measurements of /^^ have been made by several groups: CLEO-c [31 HJ |5], Belle 
[B] and BaBar [7]. Rosner and Stone combined the above measurements and report 
the experimental world average to be f^^f = (260 ± 5.4) MeV 0. Within the SM, 
has been predicted using several methods [HI El [101 ttH [El [I3l E]. While most 
calculations give values lower than the Jd^ measurement, the errors on predicted 
values are too large in most cases to claim any disagreement with experiment. The 
largest discrepancy (2.0 standard deviations) is with an unquenched lattice QCD 
(LQCD) calculation: f^^^^ = (248 ± 2.5) MeV [8]. There are several theoretical 
scenarios in which non SM particles may modify the leptonic decay rates of the Dg 
meson. Akeroyd and Chen pointed out that leptonic decay widths are modified in 
two-Higgs-doublet models (2HDM) [15]. Measurements of fn^ with an accuracy that 
matches the precision of theoretical calculations are thus necessary in order to discover 
or constrain effects of NP. 

In these proceedings we present preliminary results of absolute branching fraction 
measurements of — )■ /i^z/^ and — > T~^iyr decay^ which supersede the previous 
measurement of B{Df — )■ ^'^v^) by Belle [6] performed on a smaller data sample. 
This analysis is based on a data sample of 913 fb~^ recorded near ^/s = 10.68 GeV 
by the Belle detector at the KEKB asymmetric-energy collider [IB] . 

* Charge conjugation is assumed throughout this note unless stated otherwise. 
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2 Belle detector 



The Belle detector is a large-solid-angle magnetic spectrometer that consists of a sili- 
con vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scin- 
tillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of 
CsI(Tl) crystals located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron fiux-return located outside of the coil is instrumented to 
detect mesons and to identify muons (KLM). The detector is described in detail 
elsewhere [I7j. Two inner detector configurations were used. A 2.0 cm beampipe and 
a 3-layer silicon vertex detector was used for the first sample of 156 fb^^, while a 
1.5 cm beampipe, a 4-layer silicon detector and a small-cell inner drift chamber were 
used to record the remaining 757 fb~^ of data. 

Tracks are detected with the CDC and the SVD. They are required to have an 
impact parameter with respect to the interaction point of less than 0.5 cm in the 
radial direction and less than 1.5 cm in the beam direction. A likelihood ratio for a 
given track to be a kaon or pion, C(^k,tt), is obtained by utilizing specific ionization 
energy loss measurements in the CDC, light yield measurements from the ACC, and 
time-of-fiight information from the TOF. For electron identification we use position, 
cluster energy, shower shape in the ECL, combined with track momentum and dE/dx 
measurements in the CDC and hits in the ACC. For muon identification, we extrap- 
olate the CDC track to the KLM and compare the measured range and transverse 
deviation in the KLM with the expected values. Photons are required to have energies 
in the laboratory frame of at least 50 - 100 MeV, depending on the detecting part of 
the ECL. Neutral pion candidates are reconstructed using photon pairs with invariant 
mass between 120 and 150 MeVQ Neutral kaon candidates are reconstructed using 
charged pion pairs with invariant mass within ±20 MeV of the nominal mass. 

3 Overview of the method 

The method of absolute branching fraction measurement of — ^~T'l decays is 
similar to the one previously used by the Belle collaboration [6] and more recently by 
the Babar collaboration In this method the e^e~ — )■ cc events which contain D~ 
mesons produced through the following reactions: 

e+e- ^ cc Aag^Xfr^g^r , D*- ^ Z};7, (2) 

are fully reconstructed in two steps. In these events one of the two charm quarks 
hadronizes into a D*~ meson while the other quark hadronizes into a charm hadron 
denoted as Dtag (tagging charm hadron). We reconstruct the tagging charm hadron 

^We use natural units throughout this note. 
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as D^^ The strangeness of the event is conserved by requiring 

additional kaon, denoted as which can be either or Kg. Since Belle collected 
data at energies well above DlJ^KDl threshold additional particles can be produced 
in the process of hadronization. These particles are denoted as Xf^ag and can be: 
even number of kaons and or any number of pions or photons. In this measurement 
only pions are considered when reconstructing the fragmentation system X^ag. We 
require mesons to be produced in a D*^ -D77 decay which provides powerful 
kinematic constraint {D* mass, or mass difference between D* and Dg) that improves 
the resolution of the missing mass (defined below) and suppresses the combinatorial 
background. 

In the first step of the measurement no requirements are placed on the daughters 
of the signal Dj meson in order to obtain a fully inclusive sample of Dj events which 
is used for normalization in the calculation of the branching fractions. The number 
of inclusively reconstructed Dg mesons is extracted from the distribution of events in 
the missing mass, Mmiss(Dtag-^Xfi.ag7), recoiling against the Dtag-^Xfi.ag7 system: 

M'miss(^tag^Xfrag7) = \/ Pmiss (£'tag^Xfrag7) ^ , (3) 

where Pmiss is the missing momentum in the event: 

Pmiss(^tag^Xfrag7) = Pe+ + Pe- " PDtag ^ PK ~ PX^,,^ - P^- (4) 

Here, Pe+ and Pe- are the momenta of the colliding positron and electron beams, 
respectively, and the p_Dtag) Pk, Pxi,^ , and are the measured momenta of the 
reconstructed -Dtag, kaon, fragmentation system and the photon from Z?* — )■ Dgj 
decay, respectively. Correctly reconstructed events given in the Eq. [2] produce a peak 
in the Mmiss(-Dtag-f^Xfrag7) at nominal Dg meson mass. 

In the second step we search for the purely leptonic — )■ /x"*"//^ and — 
T^Uj. decays within the inclusive sample by requiring an additional charged track 
identified as an electron, muon or charged pion to be present in the rest of the event. 
In case of — )■ r^Ur decays the electron, muon or charged pion track identifies the 
subsequent r+ decay to e'^uj^r, fJ'^i^fJ^T or vt+Ft-. Hadronic decays, — > K^K^ and 
rjTT^ , are reconstructed partially by only requiring that an additional charged kaon 
or pion is present in the rest of the event while no requirements are set upon the 
neutral hadrons {K^ or 77) in order to increase the reconstruction efficiency. In case 
of —7- K^K^n^ all three charged tracks are required to be present in the rest of 
the event. 

■tin events where is reconstructed as tagging charm hadron additional p is reconstructed in 
order to conserve the total baryon number. 
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B [%] 




B [%] 



A+ ^ 



B l%] 



3.9 
13.9 
8.1 
4.2 
2.9 
5.4 



9.4 
6.1 
1.5 
6.9 
3.1 
1.0 



pK TT^ 



5.0 
3.4 
1.1 
1.1 
3.6 
2.6 



.0 




Sum 



38.4 



Sum 



28.0 



Sum 



16.8 



Table 1: Summary of i^tag = D^, and decay modes used in this measurement. 
The branching fractions are taken from [1]. 

4 Inclusive Ds reconstruction 

The reconstruction of the inclusive Ds sample starts with the reconstruction of the 
tagging charmed hadrons, -Dtag- In order to increase the reconstruction efficiency of 
studied events, the Z^tag is reconstructed in as many decay modes as possible, while 
keeping the purity of the sample at reasonable level. The ground state Dtag 
D~^, A+) hadrons are reconstructed in total in 18 hadronic decay modes (see Table 
[T|. Only modes with up to one 7r° in the final state are used in order to avoid large 
backgrounds. 

The center of mass momentum of the Z^tag candidates (p*) is required to be greater 
than 2.3 (or 2.5 for less clean Dtag modes) GeV/c in order to reduce the background 
and to remove charmed hadrons originating from B decays. The decay products of the 
Dtag candidates are fitted to a common vertex and candidates with poor fit quality 
are discarded {x^/n.d.f < 20). The purity of Dtag sample at this stage is rather 
low - around 17% in the signal region defined as ±3a interval around the nominal 
Dtag mass, where a is the Dtag decay mode dependent invariant mass resolution. In 
order to further clean up the Dtag sample we train an NeuroBayes \i9\ neural network 
using small sample of real data (around 1% of the total sample). Network combines 
information from the following input variables into one single scalar output variable: 
the distance between the decay and the production vertices of Dtag candidates in r — 
plane, the /n.d.f of the vertex fit of Dtag candidates, the cosine of the angle between 
the Dtag momentum vector and the vector joining its decay and production vertices 
in r — plane, the cosine of the angle between the momentum of one of the Dtag 
daughters momentum vector in the Dtag rest frame and the Dtag momentum in the 
laboratory frame (only for two-body Dtag decays), particle identification likelihood 
ratios and for the Dtag decay modes including 7r° the minimal energy of the two 
photons. To obtain the signal and background distributions of variables entering 
the NeuroBayes network a statistical tool to unfold the data distributions (sPlot) is 
applied [20] . The cut on the network output variable is optimized for each Dtag mode 
individually by maximizing S/ \/S + B, where S (B) refers to the signal (background) 
yield in the signal window of Dtag invariant mass determined by performing a fit to the 
Dtag invariant mass distribution. After the optimization the purity of the correctly 
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reconstructed Dtag candidates increases from 17% to 42% while loosing only around 
16% of signal £)tag candidates. We keep only Dtag candidates from signal region in 
-Dtag invariant mass in rest of the analysis. 

Once the ground state -Dtag hadrons have been reconstructed, and mesons 
originating from D* decays are identified by reconstructing the decays D*~^ D^n'^, 
D^n^ and D*^ — > £)°7r°, D'^'y. The motivation for this reconstruction is to clean 
up the subsequent KXfj.g^^'y reconstruction; by absorbing one more particle to the 
tagging charm hadron the subsequent combinatoric background can be reduced. In 
addition, by reconstructing D*~^ — )■ D^n^ decays we can determine of the charm 
quantum number of decays including neutral kaons. The slow pions from D* 
decay are refitted to a D production vertex in order to improve the resolution of the 
mass difference, AM — M[Dt:) — M{D). The photon energy in the laboratory frame 
is required to be larger than 50 MeV for 7r° — )> 77 and 175 MeV for D°7 
decays. In the latter decays the 7 candidate is combined with any other photon in a 
event, and if there exist a combination of two photons with invariant mass within 10 
MeV/ around the nominal 7r° mass and their energy asymmetry ( {E^-^ — E^^ ) / {E^-^ + 
Ej^)) is smaller than 0.5 the D*^ candidate is rejected. For all D* decays, the mass 
difference m{D*) — m{D) is required to be within 3(j of the corresponding nominal 
mass difference. 

A K candidate is selected to be either or K'g candidates and does not overlap 
with the L>tag candidate. 

Prom the remaining tracks and 7r° candidates in the event that do not overlap 
with the Dtg,gK candidate we form Xf^ag candidates. Only modes with up to three 
pions and up to one are used in order to avoid large combinatoric background. In 
addition pions must have momentum larger than 100 MeV/c in the laboratory frame. 
At this stage no requirement is applied to the total charge of the X^^ag system. 

The i^tag, Xfrag and K candidates are combined to form a DtagKXhag combinations 
and we keep only those with total charge ±1. The charm and strange quark content 
of the Dtg^gKXfj-g^g system is required to be consistent with that recoiling from a D*: 
if Dtag is reconstructed in ffavor specific decay mode and the primary kaon candidate 
is charged it is required that the kaon charge and the charm quantum number of 
Dtag are opposite to the D* charge; if primary kaon candidate is a Kg the charm 
quantum number of -Dtag is required to be opposite to the D* charge; and if Dtag 
is reconstructed in a self-conjugated decay mode the charge of the primary kaon is 
required to be opposite to the D* charge. A kinematic fit to Dtg^gRXf^ag candidate is 
performed in which the particles are required to originate from a common point inside 
the IP region, and the Dtag mass is constrained to the nominal value. We select only 
one DtagKXfrag candidate in an event which has missing mass, Miniss(I?tag-f^Xfrag) = 



^J\Pe+ + Pe- —PDtag ^ Pk ^PXfragP' closcst to the nominal D* mass and within 

2.00 GeV/c^ < Miniss(Aag-f^-^frag) < 2.25 GeV/c^ interval (corresponding to around 
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±3(7 interval). 

Finally, a photon candidate is identified which is consistent with the decay 
Dg'y and does not overlap with the -Dtag-ft'^frag system. We require that the energy 
of the photon candidate is larger then 120 MeV in the laboratory frame and that 
the cosine of the angle between the direction of Dtag hadron and the direction of the 
photon candidate is negative, since the signal photon should be in opposite hemisphere 
of the event with respect to the -Dtag- We perform similar kinematic fit with the 
signal photon included and with the missing mass recoiling against the Dtag-^-^frag 
constrained to the nominal D* mass. All Z)tag-^-^frag7 candidates are required to have 
p;,i,,(Dtagi^Xfrag7) > 2.8 GcV/c and M^i3s(Aagi^Xfrag7) > 1-83 GeV/c2 (see Eqs. Q 
and|3]). After the final selections, there are in average 2.1 Dtag-f^-^fragT candidates per 
event which are solely due to multiple 7 candidates. We select a best Dtag-f^-^frag7 
candidate to be the one with the highest NeuroBayes network output which is trained 
to separate signal photons from photons produced in other decays. A relative gain 
of around 23% in absolute reconstruction efficiency is obtained by applying the best 
DtagKXij-i.g'y candidate selection instead of completely random selection. Figure [l] 
shows the distribution of Miniss(Dtag-f^-^frag7) for each X^ag mode separately. 

4.1 Inclusive Ds yield extraction 

The yield of inclusively reconstructed Dg mesons is determined by performing a fit to 
the missing mass Mmiss(-Dtag-f^-^frag7) distribution for each X^ag mode individually. 
The components of the fit are divided into six categories: signal, mis-reconstructed 
signal event {K candidate or one of the pions forming Xf^ag system candidate orig- 
inate from Ds decay), D*^ — )■ D^j (signal 7 candidate originates from D*^ decay), 
D*g^ -D(s)7r° (one of the photons from vr" decay from D*^^ — >• L'(s)7r° decay is taken 
to be signal 7 candidate), wrong 7 (the energy deposited in the ECL was produced 
by unmatched charged track or by beam induced interactions), and 7 from 7r° (signal 
photon originates from a 7r° decay which does not originate from D^^^ decay). Each 
of the above six components is parameterized with a non-parametric histogram prob- 
ability density function (PDF), 7/(Mmiss(-Dtag-^-^frag7)), taken from a large sample 
of Monte Carlo (MC) simulated events. The total PDF for given Xf^ag mode is thus 
given by: 

^^*-(M^is.(Aagi^Xfrag7)) = iVsig^Jr^(M^iss(Aagi^Xfrag7)) ® e^(^^cal) 

+ E iVSk^g7^f-b4(M^i,3(Aagi^Xf,ag7)), (5) 
i=l 

where N represents the yield of each component and first (second) term in the equa- 
tion describes contribution of signal (sum of five background components). Histogram 
PDF of signal, 'Hsig'''^(Mmiss(Aag-^Xfi.ag7)), is additionally numerically convolved 
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with a Gaussian function, ^((Xcai) - centered at zero and with width acai, which takes 
into account possible differences between Miniss(-Dtag-^-^frag7) resolutions obtained on 
real data and simulated samples. 

The real data resolution in Mmiss(Dtag-^-^frag7) for signal candidates is calibrated 
using the mass difference between D* and Dg, AM = Mo* — Md^, for exclusively 
reconstructed D* — )■ Dg'y decays, where Dg decays to (pn and — )■ K^K~. In 
the exclusive reconstruction of mesons the same requirements are used for the 
signal photon candidate as in the inclusive reconstruction. The dominant contribution 
to the resolution of AM as well as of the Mjniss(-Dtag-^-^frag7) is the signal photon 
energy resolution. In former case the smearing of the Dg momentum cancels almost 
completely in the mass difference while in the latter case the impact of experimental 
smearing of Pmiss(-Dtag-ft'-^frag) on Mmiss(-Dtag-^-^frag7) is minimized by performing the 
mass constrained vertex of -Dtag-f^-'^frag candidates to the nominal mass. According 
to simulated sample the resolutions of AM and Mmiss(-Dtag-^-^frag7) are same to 
within few percent which justifies calibration of Miniss(-Dtag-^-^frag7) resolution by 
comparing AM resolutions of exclusively reconstructed D* — )■ Dg--^ decays obtained 
on real data and MC. 

In order to calibrate Miniss(-C^tag-^-'^frag7) resolution a fit to the real data AM 
distribution is performed which is described as the sum of signal and background 
components: 

J-<'^^i(AM) = iVsig?^3ig(AM) ® e;((Teai) + iVbkg(l + CiAM + caAM^). (6) 

The signal contribution, 'Hsig(AM) ® Q^Ccai), is parametrized using the histogram 
PDF obtained on simulated sample, 'Hsig(AM), and is numerically convolved with a 
Gaussian function, ^((Xcai). The background is parametrized as the 2nd order polyno- 
mial. Best agreement between real data and simulated AM distributions are obtained 
when deal = 2.0 ± 0.2 MeV. 

The Mmiss(-Dtag-ft'^frag7) distributions for each Xfrag mode of real data inclusive 
Dg sample with fit results superimposed (see Eq. |5]) are shown in Fig. [l] Total 
inclusive Dg yield on real data sample corresponding to 913 fb^^ is found to be 
94400 ± 1300, where uncertainty is statistical only. The PDFs describes well the 
observed data distributions - the normalized values of the fits are between 1.06 
and 1.32. Fit residuals exhibit no structures. To estimate the systematic uncertainty 
on the inclusive Dg yield the fits are repeated by taking cXcai = 1-8 or 2.2 MeV. 

The inclusive Dg yield, including systematics uncertainties, is found to be 

N'^''^ = 94400 ± 1300(stat.) ± 1400(syst.). (7) 

We keep only inclusive Dg candidates within the signal window of the Miniss(-Dtag-f^^frag7), 
defined as 1.95 < Mmiss(L'tagi^-^frag7) < 1.99 GeV, in rest of the analysis, except in 
case of exclusive reconstruction of Dg — )■ K~K~^7t~^ decays. 
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Figure 1: The Miniss(-Dtag-^-^frag7) distributions for all seven X^ag modes with fit 
results superimposed (solid blue line). Signal contribution is shown with the solid 
green line, while the full histograms (in different shades of gray) show the contri- 
butions of different types of background described in more details in text. The two 
dashed vertical lines indicate the signal region. 
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5 Reconstruction of Dg f decays and extraction 
of absolute B{Ds f) 

After reconstructing the inclusive sample of Ds mesons we proceed with the recon- 
struction of Ds mesons decaying to: 

• Df ^ K-K+TT+, 

• D+ K^K+, 

• ^ r]7i+, 

• D+ ^ /i+z/^, 

In the following subsections we briefly describe the reconstruction procedure and sig- 
nal yield extraction for all five studied decay modes, which then enters the calculation 
of absolute B{Ds — )■ /) as: 

B(D. /) '^f^- - 



N^:- fu.s-e{Ds^ f\md. D,) ' 



Above, N^]^^ is the number of inclusively reconstructed Dg mesons, N{Ds — J- /) is 
the number of reconstructed Ds ^ f decays within the inclusive Ds sample, and 
e{Ds /I inc. Ds) is the efficiency of reconstructing Dg f decay within the in- 
clusive Ds sample. We observed on simulated samples that the efficiency of inclusive 
reconstruction of Ds mesons depends on the Ds decay mode and therefore the inclu- 
sively reconstructed Ds sample does not represent truly inclusive sample of Ds mesons. 
The efficiency drops with increasing multiplicity of final state particles, e.g. Nch + N^^o 
(where N^h represents the number of charged particles and A^^o the number of neutral 
pions), produced in Ds decays. In order to take this effect into account a ratio of 
efficiencies to reconstruct Ds meson inclusively if it decayed to final state /, £ds-s>/) 
and the average efficiency of inclusive Ds reconstruction, e'^^ = J2i B{Ds — )■ i)e^^^^^i 
is included into the denominator of Eq. |8j /bias = ^^ol-^f I^^dI- Ratio /bias is taken 
from the simulated sample including all known Ds decay modes. 

Possible differences between simulated and real data samples in terms of Ds decay 
modes and their i3s used in our simulation are estimated by studying the distribu- 
tions of number of particles, A'ch + A'ttO, produced in Dg decays. This distribution 
on real data is obtained in the following way: first for each inclusive Ds candidate 
the number of remaining charged tracks, Nl^°, and vr" candidates, N^l'^°, not asso- 
ciated to Dtag-f^A'fi.ag7 candidate are counted; second fits to the M^iss{Dta,gK Xf^g^^'y) 
distributions are performed in bins of A^h'^" + A^o"^" and the inclusive Ds signal yield 
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is recorded. The distribution of N^^^° + N^l^° is roughly proportional to the true 
distribution of N^h + N^^o, but with a considerable amount of convolutioiij^ The un- 
folded Nch + N^o distributions are obtained from N^^^° + A^^o™ distribution using 
the singular value decomposition algorithm. From the distributions of A'ch + N^^o 
obtained on real data and simulated samples and the dependence of inclusive Dg 
reconstruction efficiency obtained on simulated sample we estimate the the ratio be- 
tween real data and simulated average inclusive Dg reconstruction efficiencies to be: 
Idata/^Ss Imc = 0.9768 ± 0.0134. The ratio is consistent with unity within the 
uncertainty. Nevertheless the inverse of the /bias factor is corrected by the above ra- 
tio and its error is taken as a source of systematic uncertainty of measured absolute 
branching fractions. 

5.1 Df K-K+7T+ 

The reconstruction of — )■ K^K^tt^ decays is performed by requiring exactly 
three charged tracks in the rest of the event with a net charge equal to the charge 
of inclusively reconstructed Ds candidate. The track with opposite charge to the 
inclusive D'^ candidate is selected to be K~ candidate while the two same-sign tracks 
are identified as or vr"*" based on their likelihood ratios, Ck,tt- 

The exclusively reconstructed Df — )■ K'K'^tt^ candidates are identified as a 
peak at the nominal mass of the D*^ in the invariant mass distribution of K^K^tt^'j 
combination, M{K^K^'k^^). Here the 7 stands for the signal photon candidate 
used to reconstruct -Dtag-^-^frag7 candidate. The M{K^ K~Ti^'-)) is chosen over 
M{K^K~Ti^)^ because both sides (inclusive and exclusive) have to be correctly re- 
constructed to produce peaks in M{K^ K'n^--^) and Mmiss(-Dtag-f^-^frag7)- Correctly 
reconstructed — >■ K~K^'k^ events will namely peak in M{K^K~tt^) even if the 
inclusive reconstruction of Ds candidates failed, e.g. the photon candidate which 
does not originate from — )■ decay is picked up. Since M{K~ K^n^'-^) and 
Mmiss(Dtag-^-^frag7) are correlated due to their common input - 7 4- momentum - no 
additional selection is applied to the Miniss(Dtag-^-^frag7)- 

The signal yield of exclusively reconstructed Dg — ?■ KKn decays is extracted by 
fitting the M{K^K^7T^'~f) distribution. The components of the fit are divided into 
three categories: correctly reconstructed D* Dg'-f — )■ KKti^ decays (signal); mis- 
reconstructed Dl — )• DgTT^ — )■ K^K~Ti^ decays, where one of the photons from 
the 7r° decay is not reconstructed {D* — )■ Dgir^ — )• KKtt'j^); and random combina- 
tions of charged tracks and photons (combinatorial background). The M[K~^K~7r~^^) 
distribution is parameterized as: 

^E.g. Dg daughter particle might not be reconstructed or a fake charged track or tt^ candidate 
is counted. 
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+iVcomb ■ (1 + ci ■ M{K+K-n+^) + 

C2 ■ M{K+K~n+-fy + C3 ■ M{K+K-TT+-ff), (9) 

where the signal and D* — > Dgir^ background are parametrized using the non- 
parametric histogram PDF obtained on a simulated sample (7{), and the combinato- 
rial background is parametrized as a third order polynomial. The signal distribution is 
convolved with a Gaussian function, ^(cr^^f^), in order to take into account the differ- 
ences between the resolutions of M^K^K'n'^'y) on real data and simulated samples. 



The (j^^f^ is estimated using the same procedure as described in section 4.1 The only 
difference is that the resolution on M{KK'K'y) is calibrated instead of the and Dg 
mass difference. We determine a'^^ = 3.2 ± 0.2 MeV. Free parameters of the fit are 
the normalization parameters, Ni, and combinatorial background shape parameters, 

Ci. 

The M{K~^K~7r^'-f) distribution of exclusively reconstructed Dg — )■ KKn decays 
within the inclusive Dg sample obtained on real data sample is shown in Fig. |2] with 
fit results superimposed. The number of correctly reconstructed Dg — )■ KKn decays 
is found to be 

N{Dg KKtt) = 4094 ± 123, (10) 
where the error is statistical only. 



5.2 D+ K^K+ 

The — 7- K^K^ decays are reconstructed partially by requiring only one additional 
charged track to be present in rest of the event. The charged track is required to 
be consistent with kaon hypothesis and has charge equal to that of the inclusively 
reconstructed Dg candidate. The neutral kaon is not reconstructed at all and it is 
identified as a peak at the nominal mass of neutral kaon squared in the missing mass 
squared distribution: 

where the missing 4-momentum is given by 

Pmiss(£'tagi^^frag7^) = Pe+ + Pe" " PDtag - PK - PXf,,g -p^f-pK- 

An explicit reconstruction of K'^ meson (via the experimentally accessible Kg 
TT+TT^ decays) would lead to a significant signal loss (2/3 of the signal would be lost, 
without even taking into account the Kg reconstruction efficiency). The signal peak 
of Df — )■ K^K~^ in M^^Dts.gKXfj-g^g'yK) distribution is in addition used to study 
the differences between M^jgg resolutions in simulated and real data samples, which 
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Figure 2: The M{K'^K~7r^'~f) distribution of exclusively reconstructed Dg — ?■ KKn 
decays within the inclusive Dg sample with fit results superimposed (solid blue line). 
Solid green line shows the signal contribution, while the full histograms show the 
contribution of combinatorial background (in red) and Dg — )■ KKtt decays originating 
from D* — 7- DgTc^ decay (in blue). 

is important in the extraction of yield of — )■ /x"*"//^ decays. Since the flavor of the 
neutral kaon is not determined, the doubly Cabibbo suppressed decays, — )■ K^K^ , 
also contribute to the peak in M^jgg(Dtag-f^-^frag7-^)- Their relative contribution can 
naively be estimated to be equal to tan^ 6'c7 ~ 0.29% {6c being the Cabibbo mixing 
angle), which is order of magnitude below the expected statistical uncertainty and 
can thus be safely neglected. 

The signal yield of partially reconstructed — )■ K^K~^ decays is extracted by 
performing a fit to the M^igg(L'tag-f^-^frag7-^) distribution. The candidates are divided 
into six categories: candidate originates from the K^K^ decay and the 

inclusive Dg candidate is correctly reconstructed (signal); candidate originates 
from the — )■ rjK^ decay and the inclusive Dg candidate is correctly reconstructed 
{Dg — 7- rjK decays); candidate originates from the — )► -k^K^ decay and the 
inclusive Dg candidate is correctly reconstructed {Dg — )• tx^K decays); candidate 
is mis-reconstructed pion originating from the D'^ — )■ r^vr^ decay and the inclusive 
Dg candidate is correctly reconstructed {D^ — )■ t^tt^); candidate originates from 
the Dg K*^K^ K^Tx^K^ decay and the inclusive Dg candidate is correctly 
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reconstructed {Ds — t- K*^K^ decays); and all other candidates (combinatorial back- 
ground). The M^jgg(Dtag-^-^frag7-^) distribution is parameterized as: 

:F{MI,,,) = iV,ig.^^(M^j,„m^o,3.af^) 

i=l 

3 3 
i=l i=l 

+iVcomb ■ (1 + ClM^.^, + C2M2.^, + CsM^;^, + C4M;^i,J, (11) 

where the signal, and the two peaking backgrounds {Ds — )■ rjK and Dg — )■ Tr^i^) 
are parametrized using the sum of 3 Gaussian functions. All the parameters of the 
latter peaking backgrounds are fixed to values determined on the simulated sample. 
In case of the signal all the shape parameters are fixed to the values determined on 
simulated sample, except the mean (wt-I-o) and the common resolution scaling factor 
of the core and the second Gaussian (si = S2 = s). The width of the third Gaussian 
(describing the outliers) is fixed. The rjir and K*K^ candidates are parametrized 
with a sum of Gaussian and a Bifurcated Gaussian function, where all parameters are 
fixed to the values determined on simulated sample. The combinatorial background 
is parametrized with a polynomial of the 4*^^ order, where the coefficients Cj are de- 
termined with the fit to the M^jgg(Dtag-^-^frag7-^) distribution for candidates in the 
^miss(-Dtag-f^-^frag7) < 1-95 GcV sidcbaud region. Yields of all event categories are 
free parameters of the fit, except for the Ds — t- rjK and Ds — )■ r/vr, which are con- 
strained to expected values based on their known Bs and MC determined efficiencies. 

The M^igg(Aag-f^-^frag7-f^) distribution of partially reconstructed -> K°K~^ 
decays within the inclusive Ds sample obtained on real data sample is shown in Fig. 
[3] with fit results superimposed. The number of correctly reconstructed D'^ — )■ K^K~^ 
decays is found to be 

N{Ds ^ K°ir) = 1943 ± 82, (12) 
where the error is statistical only. 



5.3 Df rj7r+ 

As in case of Df — )■ K^K^ decays we perform partial reconstruction of D"^ — ?► rjn^ 
decays as well. We require only one charged track consistent with the pion hypothesis 
to be present in rest of the event. We do not perform any explicit reconstruction of r] 
mesons which are identified as a peak at the nominal mass squared of rj in the missing 
mass squared distribution: 

^miss(Aag^^frag77r) = (^tag^^frag77r) , 
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Figure 3: The M^^^^^Dts.gKXij-ag'yK) distribution of partially reconstructed D]. 
K^K^ decays within the inclusive Dg sample with fit results superimposed (solid 
blue line). Solid green line shows the signal contribution, while the full histograms 
show the contribution of combinatorial background (in red), — > tt^K^ (dark gray), 
— )■ TjK^ (light gray) or other decays (blue). 



where the missing 4-momentum is given by 

Pmiss(^tag^^frag77r) = Pe+ + Pe- " P^ag - PK - PX^.g " P7 " P^' 

An explicit reconstruction of 7] meson would lead to a significant signal loss. 

In this sample of inclusive Ds candidate plus additional charged pion there is a 
significant contribution from — )■ rv decays, when r lepton decays hadronically 
to charged pion and neutrino. These events are suppressed by requiring that the 
extra neutral energy in the ECL (-Eecl) to be larger than 1.0 GeV, where the -Eecl 
represents a sum over all energy deposits in the ECL which are not associated to the 
tracks or neutrals used in inclusive reconstruction of Ds candidates nor the charged 
pion candidate [21j. The Dg ^ tu ^ nuu decays namely peak in Eecl at 0, while 
Dg — )■ rjTr decays deposit significant amount of energy in the ECL via the rj decay 



products (see section 5.5 for more details). 

The signal yield of partially reconstructed — ?■ rjir^ decays is extracted by 
performing a fit to the M^igg(Dtag-f^-^frag77r) distribution. The candidates are divided 
into five categories: 7r+ candidate originates from the Df — )■ r^vr"'' decay and the 
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inclusive Ds candidate is correctly reconstructed (signal); 7r+ candidate originates 
from the — )■ K^ix^ decay and the inclusive Dg candidate is correctly reconstructed 
[Dg —7- K^Tt decays); vr^ candidate is mis- reconstructed kaon originating from the 
— )■ K^K~^ decay and the inclusive Dg candidate is correctly reconstructed {Dg — )■ 
K^K decays); 7t~^ candidate originates from the Dg — )■ p^K^ — )■ tt^h'K^ decay and 
the inclusive Dg candidate is correctly reconstructed {Dg — )• p^K^ decays); and all 
other candidates (combinatorial background). The M^igg(Dtag-^-^frag77r) distribution 
is parameterized as: 

i=l i=l 

+iVcomb ■ (1 + ClM^j^, + + csM^,, + c,M^,J, (13) 

where the signal and the peaking background from Dg — > K^tt decays are parametrized 
using the sum of three Gaussian functions. All the parameters of the latter peaking 
background are fixed to values determined on simulated sample. In case of the signal 
all the shape parameters are fixed, except the mean (m^). The common resolution 
scaling factor of the core and the second Gaussian [si = S2 = s) is constrained to the 
value obtained from the fit to the Dg — )■ K^K sample. The width of the third Gaus- 
sian (describing the outliers) is fixed. The K^K and pK candidates are parametrized 
with a Bifurcated Gaussian function, where all parameters are fixed to the values de- 
termined on simulated sample. The combinatorial background is parametrized with a 
polynomial of the 4*'' order, where the coefficients Cj are determined with the fit to the 
M^igg(L'tag-f^-^frag77r) distribution for candidates in the Mmiss(-Dtag-f^^frag7) < 1-95 
GeV sideband region. Yields of all event categories are free parameters of the fit, 
except for the K^n and K'^K, which are constrained to expected values based on 
their known Bs and efficiencies determined on simulated sample. 

The M^jgg(Z}tag-^-^frag77r) distribution of exclusively reconstructed D'^ — )• 7771+ 
decays within the inclusive Dg sample obtained on real data sample is shown in 
Fig. [4] with fit results superimposed. The number of correctly reconstructed Df — 
T^TT"'" decays is found to be 

N{Dg riTi) = 773 ± 58, (14) 
where the error is statistical only. 

5.4 D+ ^i+i^^ 

The Df — )• /i^z/^ decays are reconstructed by requiring additional charged track 
consistent with muon hypothesis to be present in rest of the event. Single missing 
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Figure 4: The Mj^jgg(Dtag-^-'^frag77r) distribution of partially reconstructed — )• rjir^ 
decays within the inclusive Dg sample with fit results superimposed (solid blue line). 
Solid green line shows the signal contribution, while the full histograms show the 
contribution of combinatorial background (in red), — )■ K'^n'^ (dark gray), or 
other decays (blue). 



neutrino is then identified as a peak at zero in the missing mass squared distribution: 

A^miss(Aag^Xfrag7/^) = Pmiss (^tagi^^fragT/^) , 

where the missing 4-momentum is given by 

Pmiss(^tag^^frag7/^) = Pe+ + Pe" " P^ag - PK - PXfr.g " P7 " 

The signal yield is extracted by performing a fit to the M^igg(Dtag-^-^frag7yu) dis- 
tribution. The candidates are grouped into five categories: /x"*" candidate originates 
from the — )■ /i"^z/^ decay and the inclusive Dg candidate is correctly reconstructed 
(signal); fi~^ candidate is mis- reconstructed kaon originating from the — )■ K^K~^ 
decay and the inclusive Dg candidate is correctly reconstructed {Dg — )■ K^K decays); 
/i"^ candidate is mis-reconstructed pion originating from the — > rj-K^ decay and 
the inclusive Dg candidate is correctly reconstructed {Dg rjir decays); fi~^ candi- 
date originates from the Dg rur ^i'tVi_i decay and the inclusive Dg candidate 
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is correctly reconstructed {Dg — ?■ rur); and all other candidates (combinatorial back- 
ground). The M^igg(Z)tag-^-^frag7/^) distributions is parameterized as: 

3 2 

1=1 i=l 

+NrExp(tr) + A^comb-Exp(tcomb), (15) 

where the signal and the peaking background Dg — t- rj-K are parametrized using the 
sum of 3 and 2 Gaussian functions, respectively. All the parameters of the latter 
peaking background are fixed to values determined on MC sample. In case of the 
signal all the shape parameters are fixed, except the mean (mg). The common reso- 
lution scaling factor of the core and the second Gaussian (si = S2 = s) is constrained 
(not fixed) to the value obtained on the Ds K^K and Ds i]7c samples. The 
width of the third Gaussian (describing the outliers) is fixed. The K^K candidates 
are parametrized with a Bifurcated Gaussian function, where all parameters are fixed 
to the values determined on MC sample. The Dg — )• tUt- and combinatorial back- 
grounds are parametrized with a exponential function, where the shape parameter of 
the combinatorial background is determined with the fit to the M^jgg(L'tag-ft'-^frag7Ai) 
distribution for candidates in the Mniiss(-Dtag-ft'^frag7) < 1-95 GeV sideband region 
and it is fixed to MC determined value for riyr decays. Yields of all event categories 
are free parameters of the fit, except for the rj7r and K^K, which are constrained to 
expected values based on their known Bs and MC determined efficiencies. An im- 
portant peaking background could arise from the — )■ tt+tt^ decays, which would 
peak near the signal fJ'^'^n decays (at m^o = 0.018 GeV^). These decays were 

not observed so far and based on the upper limit of B{Df — )■ tt+tt^) < 3.4 x 10~^ at 
90% C.L. set by CLEO-c [22j their contribution is estimated to be less than 1 event. 
Their contribution is therefore negligible and is not taken into account. 

The distribution of Mmiss(-Dtag-f^-^frag7yu) with superimposed fit result is shown in 
FigjS] The number of reconstructed Dg — )■ /iz/ decays is found to be 

N{Dj /i+z/^) = 489 ± 26, (16) 

where the error is statistical only. 

5.5 D+ T+l^r 

The reconstruction of t^Vt is performed by requiring one charged track to be 

present in rest of the event and is identified as an electron, muon or a pion (denoted 
as Dg — )• t{X)u where X = e, /x, vr) indicating subsequent decay of r lepton to ei/gi^r, 
Ixp^Pr or ttz/t-. Due to the multiple neutrinos in the final state these decays do not 
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Figure 5: The M^igg(Dtag-f^-'^frag7/^) distribution of exclusively reconstructed D"^ — )► 
yU"*"!/^ decays within the inclusive Ds sample fit results superimposed (solid blue line). 
The solid green line shows the contribution of signal, while the contributions of com- 
binatorial background and background other Ds decays are indicated by the shaded 
red and blue areas, respectively. 

peak in the missing mass squared distribution: 

^miss(^tagi^^frag7^) = (^tagi^Xfrag7X) , 

where the missing 4-momentum is given by 

Pmiss(Aag^Xfrag7^) = Pe+ + Pe" " P^tag - PK - PXf,,g -p^-pX- 

The background in — r(7r)z/ sample is much larger than in the leptonic modes 
which however can be significantly reduced by requiring the missing momentum of 
the event, Pmiss(-Dtag-^-^frag77'")) to be larger than 1.2 GeV. Background in this sample 
can be further reduced by requiring 0.0 < M^igg(Z)tag-^-'^frag77r) < 0.6 GeV^ In case 
of Dg — )■ TUr — )■ iur^r^e wc require M^^^^{DtagKXi^g,g1^) > 0.3 GeV^ in order to veto 
contribution of — )► i^Ui decays. 

^Due to the lack of phase space and the fact that Dg — > Ti^r — > t^i^t decays have only two neutrinos 
in the final state these decays populate relatively narrow region in M^-^^^{DtiigK Xf^i^g'^rr) . 
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The signal yield of Dg T{X)vr decays is extracted from the fits to the Eecl 
distributions. The candidates are divided into three categories: X originates from the 
Dg — )■ r(X)z/^ decay and the inclusive Dg candidate is correctly reconstructed (signal); 
small contribution from other r decays (cross-feed) is also considered as signal; X 
originates from a Dg decay such as semileptonic hlv decays (when X = £) or hadronic 
Dg decay modes (when X = ti) and inclusive Dg candidate is correctly reconstructed 
{Dg —7- / background); and all other candidates (combinatorial background). The 
Eecl distribution is parametrized as: 

•^(-E'eCl) = ^r(X) ■ '^r(X)(-£^ECL) 

+ X! ^Ds^f ■ 'HDs^fiE-ECh) + Ai'comb " 'Hcomb(-£'ECL), (17) 
/ 

where all components are described with non-parametric histogram PDF ("H) taken 
from simulated samples. 

The distributions of -Eecl for Dg — )■ r(e)z/, Dg — )• t{ii)v, and Dg — )■ r(7r)z/ decays 
are shown in Fig. |6] with fit results superimposed. The number of reconstructed signal 
decays are found to be: 

N{Dg T{e)ur) = 952 ± 59, (18) 
N{Dg-^T{fx)ur) = 758 ±48 (19) 
N{Dg T{n)ur) = 496 ± 35, (20) 

where the errors are statistical only. 

6 Results and Systematics 

From the extracted signal yields of studied Dg decay modes we determine their ab- 
solute branching fractions using Eq. [8j They are summarized in Table |2] and are 
found to be consistent with previous measurements performed by CLEO and Babar 
collaborations. 

Systematic errors for the measured branching fractions are associated with the 
uncertainties in the signal yields, the efficiencies, and the number of inclusively re- 
constructed Dg. 

The systematic erroi{|i] related to the normalization is assigned to be 1.95% (see 
Eq. [T]) and is common for all studied Dg ^ f decays. As described in Sec. |5] the 
/bias factor is corrected by the FJ^^'Idata/^Ss'Imc ratio and the uncertainty on the 
correction (1.37%) is assigned as systematic error referred to as tag bias and is com- 
mon for all Dg ^ f decays. The systematic errors in the Dg ^ f reconstruction 

"All systematic errors are given as relative uncertainties. 
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Figure 6: The -Eecl distribution of exclusively reconstructed Dg — )■ T{e)Vr (top left), 
Dg — > t{^)vj. (top right), and Dg — )■ t{tt)vj. (bottom) decays within the inclusive Dg 
sample with fit results superimposed. The solid green lines show the contribution 
of signal and the contribution of r cross-feeds are indicated by the full dark green 
histograms. The contributions of combinatorial background and background from 
other Dg ^ f decays is indicated by the shaded red and blue areas, respectively. 



efficiencies arise from the uncertainty in tracking efficiency (0.35% per reconstructed 
charged track in the final state of Dg ^ f decay), particle identification efficiency, 
branching fractions of r decays, and MC statistics. All of these depend on the stud- 
ied Dg ^ f decay. The systematic error from MC statistics of the histogram PDFs 
is evaluated by varying the content of each bin by its statistical uncertainty in case 
of Dg — > TUr signal yields. The shape parameters of combinatorial background are 
varied within their uncertainties and the differences with respect to the nominal fits 
are taken as the systematic uncertainty in case of Df — )• n'^i'^, Dg — )■ K^K^^ and 
Dg — )■ T^TT decays. To estimate the systematic error due to the possible signal -Eecl 
shape difference between MC and data, the ratio of data to MC for the background 
subtracted -Eecl histograms of the Dg — )■ KKtt and Dg — )■ K^K^ samples is fit- 
ted with a fourth-order polynomial and the signal Eecl PDF is modified within the 
fitted errors. Similarly we estimate the systematic error due to the possible combi- 
natorial Eecl shape between MC and data by fitting the ratio of Eecl histograms of 
Dg — )■ Tu^ candidates populating the Mfniss(-Dtag-^-^frag7) sideband region (defined as 
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Z)+ Decay Mode Signal Yield /bias ■ ^{Ds /|incl. 



B 



4094 ± 123 
1943 ± 82 
773 ± 58 



0.8578 ± 0.0056 
0.7254 ± 0.0054 
0.4577 ± 0.0064 



5.06 ±0.15 ±0.19 
2.84 ±0.12 ±0.08 
1.79 ±0.14 ±0.05 



489 ± 26 



0.9813 ±0.0175 



0.528 ±0.028 ±0.019 



5.37 ± 0.33^^:;^^ 

5.88 ± 0.37i[JJ^ 
5.96 ± QA2tlfg 
5.70±0.2ll[;:ii 



TUr (e mode) 
(/^ mode) 
TUj. (tt mode) 
TVr (Combined) 



952 ± 59 
758 ± 48 
496 ± 35 



1.0531 ±0.0135 
0.7849 ±0.0118 
0.8072 ±0.0152 



Table 2: Signal yields, tag bias corrected efficiencies and measured branching fractions 
for all five studied decay modes. The ffist uncertainty is statistical and the second 
is systematic. 



Mmiss(Aagi^^frag7) < 1-94 GcV or Mmi,,(Aag^^frag7) > 2.00 GeV) with a linear 
function and modifying the combinatorial i?ECL PDF within the fit errors. The un- 
certainties for the branching fractions of Dg decays the true Dg background categories 
that are fixed in the fits are estimated by changing the branching fractions by their 
experimental errors [1]. In case of Dg — )■ T{e)Vr and Dg — t- T{^)vr decays the largest 
contribution to the systematic uncertainty originates from Dg — t- K^lv decays and in 
case of Dg — )■ t{tt)vj. from Dg tjtt and pK decays. The r cross-feed is fixed relative 
to the signal contribution in the nominal fit to the E-^cl distributions of Dg — )■ tu 
samples. The ratios are varied within their uncertainties and the fits are repeated 
and the differences from the nominal fits are taken as the systematic uncertainties. 

The total systematic error is calculated by summing the above uncertainties in 
quadrature. The estimated systematic errors are summarized in Tables |3] and |4] for 
the hadronic and leptonic Dg decay modes, respectively. 



7 Extraction of f^^ and Conclusions 

The value of /^^ is determined from measured branching fractions of leptonic Dg 
decays. Inverting Eq. [T] yields 

8nB{Dg -> iui) 

The external inputs necessary in the extraction of jr,^ from the measured B's, are given 
in Tablejs} The |Vcs| is obtained from the very well measured = 0.97425(22) and 
iKbl = 0.04 from an average of exclusive and inclusive semileptonic B decay results 
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\%\ K^K^ \%] 




Statistical 


o.UU 


4.ZZ 


r .i3U 


INForm ?i 1 i 7 ^1 1 i on 




1 Q5 


1 Q5 


Tag bias 


1.37 


1.37 


1.37 


Tracking 


1.05 


0.35 


0.35 


Efficiency 


0.65 


0.74 


1.40 


PID 


2.57 


0.82 


1.08 


Signal PDF 


0.82 






True Ds background 




0.56 


0.65 


Total syst. 


3.81 


2.71 


3.06 


Stat. + Syst. 


4.85 


5.02 


8.10 



Table 3: Summary of systematic uncertainties for the branching fraction measurement 
of hadronic Dg decays. 



Source 



Tie p 7o 



t{^)v [%] t{'k)v 



TV r/o 



Statistical 


±5.32 


±6.18 


±6.33 


±7.04 


±3.75 


Normalization 


±1.95 


±1.95 


±1.95 


±1.95 


±1.95 


Tag bias 


±1.37 


±1.37 


±1.37 


±1.37 


±1.37 


Efficiency 


±1.78 


±1.28 


±1.51 


±1.88 


±0.84 


Tracking 


±0.35 


±0.35 


±0.35 


±0.35 


±0.35 


PID 


±1.96 


±2.03 


±1.93 


±0.88 


±1.70 


Signal PDF 




±3.46 


±1.96 


±3.43 


±2.95 


Comb. bkg. PDF 


±0.02 


±0.11 


-8.31 


±0.92 


-2.54 


PDF Stat. 




±2.16 


±2.19 


±3.05 


±1.44 


True Ds background 


±0.82 


±3.88 


±3.56 


±3.15 


±2.84 


T cross-feed 




±0.36 


±0.24 


±3.71 


±0.94 


B{t X) 




±0.22 


±0.23 


±0.64 


±0.19 


Total syst. 


±3.67 


+6.59 
-5.61 


+5.76 
-9.92 


+7.49 
-6.60 


+5.40 
-5.19 


Stat. + Syst. 


±6.46 


+9.03 
-8.35 


+8.56 
-11.8 


+ 10.3 
-9.65 


+6.57 
-6.40 



Table 4: Summary of systematic uncertainties for the branching fraction measurement 
of leptonic Dg decays. 



as discussed in Ref. [23] by using the following relation: 

1, 



IK. 



IK 



ud 
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Quantity 



Value 



1.96847(33) GeV 
1.77682(16) GeV 
0.105658367(9) GeV 
0.500(7) ps 
1.16637(1) X 10-5 GeV-2 
0.97345(22) 

Table 5: Numerical values of external parameters used in extraction of Jds- 

Ds ^ /d. [MeV] 



TV 


249.0 ± 6.6(stat.; 
261.9 ±4.9(stat.; 


) ±4.6(syst.; 
) ± 7.0(syst.; 


)±1-7(^dJ 
)±l-8(r,,J 


Combination 


255.0 ±4.2(stat.; 


)±4.7(syst.; 


)±l-8(r,5j 



Table 6: Measured values of /d^ in (iv and tv decay modes and the combination of 
the two. 

The external inputs are all very precisely measured and do not introduce additional 
uncertainties except the Ds lifetime, tds, which introduces an 0.70% relative uncer- 
tainty on fr,^. Table [6] summarizes the obtained values of fn^ using the fJ-^i^n 
and — t- r+z/,- decays. An error-weighted average of /^^ is found to be 

/^^ = (255.0 ± 4.2(stat.) ± 4.7(syst.) ± 1.8(rDj) MeV, 

where the correlation of the systematic uncertainties between the fiv and tv have 
been taken into account. This is the most precise measurement of fn^ up to date. 
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